STABILITY OF L°° SOLUTIONS FOR HYPERBOLIC SYSTEMS WITH 
COINCIDING SHOCKS AND RAREFACTIONS 



STEFANO BIANCHINI 



Abstract. We consider a hyperbolic system of conservation laws 

( u t + f{u) x = 
\ u(0, ■) = «o 

where each characteristic field is either linearly degenerate or genuinely nonlinear. Under the assumption 
of coinciding shock and rarefaction curves and the existence of a set of Riemann coordinates w, we prove 
that there exists a semigroup of solutions u(t) = Stuo, denned on initial data uq £ L°° . The semigroup 
S is continuous w.r.t. time and the initial data Uq in the Lj L oc topology. Moreover S is unique and its 
trajectories are obtained as limits of wave front tracking approximations. 



S.I.S.S.A. Ref. 65/2000/M 



(1.1) 



1. Introduction 

Consider the Cauchy problem for a strictly hyperbolic system of conservation laws 

u t + f(u) x = 



u(0, •) = u 



where u £ K" and / : Vt \— > R" is sufficiently smooth, 51 open. If the initial data it is of small 
total variation, the global existence was proved first in Jr7fl . Moreover a series of paper ^, || ||, [l4j 
establishes the uniqueness and well posedness of the Cauchy problem ( |1.1| ). However, when uq has large 
total variation or even more generally uq belongs to L°°, the solution u may not exist globally in L°° 
^9| : only for special system it is possible to consider initial data with large total variation. We recall 
some of the results available in this direction. 



1) For scalar conservation laws, the entropic solution to (1.1) generates a contractive semigroup w.r.t 
the L 1 distance, on a domain of L°° data fl20f . 

2) For general Temple class system, in [0, [|, |23| it is proved the existence and stability of the entropic 
solution for initial data with arbitrarily large but bounded total variation. 

3) If all characteristic families are genuinely nonlinear and the system is Temple class, the existence 
and stability for initial data in L°° is proved in [ pd| . 

4) For special 2x2 systems, in which one of the equation is autonomous, various results have been 
proved in [IB) , with initial data with unbounded total variation. 

An open question is if the semigroup of solutions to the systems of case 2), defined on all the initial data 
Uq with total variation arbitrary large but bounded, can be extended to data in L°°. In many model 
systems, in fact, some of the characteristic fields are linearly degenerate, so that the result in |ll| does 
not apply: for example, the traffic model considered in Q] or a 2 x 2 model for chromatography. Aim 
of this paper is to prove that, at least in the case where the flux function / is convex and shocks and 



rarefactions coincide, the solution to (1.1) can be defined for uo 6 L 



To illustrate the heart of the matter, we assume that the system (1.1) admits a system of Riemann 
coordinates w £ R", and that shock and rarefaction curves coincide in f2. Moreover we assume that 
each characteristic field is linearly degenerate or genuinely nonlinear. We do not assume that rarefaction 
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curves are straight lines. We consider a set E of the form 

E = I it G fl : w(u) G [ai, 6i], i = 1, . . . , n\. 

With L°°(R; E) we denote the space of L°° functions with values in E. The main result of this paper is 
the following: 

Theorem 1.1. There exists a unique semigroup S : [0, +oo) x L°°(R;E) i — ► L°°(M.;E) such that the 
following properties are satisfied: 

i) for all u n , u G L°°(R; E), t n , t G [0, +oo), with u n —* u in L\ oc , \t — t n \ — > as n — > +oo, 

lim S tn u n = S t u in L] ■ 

n — »+oo 



ii) each trajectory StUo is a weak entropic solution to the Cauchy problem (1.1) with uo G L°°(R;-E); 
hi) if uq is piecewise constant, then, for t sufficiently small, StUQ coincides with the function obtained 
by piecing together the solutions of the corresponding Riemann problems. 



The first two properties imply the well posedness of the Cauchy problem (1.1 ) with initial data in L°°. 
The uniqueness follows because S satisfies iii) and it is limit of wave front approximations. 

As it is shown in the last example of pi| , the semigroup S cannot be uniformly continuous: thus 
we cannot apply any compactness argument to construct the solution u(t) = StUQ. The fundamental 
problem is that, differently from [pi, the total variation of the Riemann invariants corresponding to 
linearly degenerate families does not decrease in time. 

The main idea of this paper is to study how the solution to the characteristic equation 

(1.2) x(t) = X t (u(t,x(t))), x(Q)=y, 



depends on the solution u of (IT). It will be shown that, for a fixed time r, the map y i— > x(t) depends 
Lipschitz continuously on the initial data uq, and moreover the Lipschitz constant is independent of the 
total variation of uo- Since the Riemann invariant Wi is the broad solution to 

{wi) t + Xi{u(t, x))(vii) x = 0, 

a simple argument gives the convergence of the wave front tracking approximations. Note that this also 



implies the well posedness of the ODE ( 1.2) wh en u(t, x) is an L°° solution of the system (IT). As noted 



in [Kj, for general system the solution to (1.2) does not exist. 



The paper is organized as follows. Section |^ contains the basic assumptions on the system (1.1) 



Moreover we construct the wave front approximation of the solution u(t). In Section || we analyze 
carefully the shift differential map, i.e. the evolution of a perturbation in uq in which only the position of 
the initial jumps has changed. The method we use is essentially the one in with slight modifications 
due to the fact that in our system the rarefaction curves do not need to be straight lines. The main result 
is here the explicit computation of the shift differential map. 



Section ^ is concerned with the equation for characteristics (1.2). We prove the Lipschitz dependence 



of the map y i— > x(t) w.r.t. both the initial data uo and y. Moreover we will show that the Lipschitz 
constant is independent from the total variation of uq. Finally, in Section ^, we prove Theorem [l.l[ 

2. Basic assumptions and wave front approximations 
We consider a strictly hyperbolic system of conservation laws 
(2.1) «f + /(«)x = 0, 

where / : f2 — > W 1 is a smooth vector field defined on some open set ft C R". Let A(u) = DF(u) 
be the Jacobian matrix of / and denote with Xi(u) its eigenvalues and with ri(u), l l (u) its right and 
left eigenvectors, respectively. We assume that the eigenvalues Xi can be either genuinely nonlinear or 
linearly degenerate. In the following the i-th rarefaction curve through u G £1 will be written as Ri(s)u, 
with Ri(0)u — u, while the i-th shock curve will be denoted by Si(s)u, and its speed by o~i(s,u). The 
directional derivative of a function <f>{u) in the direction of r, (it) will be denoted as 

<f>(u + hn(u)) -(f>(u) 

n • 4>(u) = hm , 

h->0 h 
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while the left and right limit of a BV function / in a point x will be written as 

f(x-)= lim f(y), f(x+)= lim J{y). 

y — >x— y — >x+ 

We assume that the rarefaction curves Ri generate a system of Riemann coordinates w(u). We recall 
that a necessary and sufficient condition for the local existence of Riemann coordinates is the Frobenius 
involutive condition: if [X,Y] denotes the Lie bracket of the vector fields X,Y, the condition is 

Tj] G span{ri,rj} for all i,j — 1, . . . ,n. 

In the following we will use indifferently the conserved coordinates u or the Riemann coordinates w. 
Fix a domain 

(2.2) E = {u G O : w(u) G [a^h], i = 1, • . . , n|. 
Since is compact, there is a constant c > such that 

(2.3) Ti • Xi(u) > c \/u £ E, if A,; is genuinly nonlinear. 



We suppose that the system (2.1) is uniformly strictly hyperbolic in f2: this means that there exists a 
constant d such that 

(2.4) X i+1 (u) - Xi(v) > d, Vu,v G E, i = 1, .. . ,n- 1. 



We also assume that in the system ( |2.1[ ) shock and rarefaction curves coincide: this implies |26J that 
either the rarefaction curve Ri(s)u is a straight line or the eigenvalue is linearly degenerate. In fact, one 
can prove that 



(2.5) ^ (s ' u) 



= rfc • \i{u))(l l (u),ri • r f (u)) + -r-j • (r< • A;(u)), 

s=0 



and for the shock curve Si(s)u we have 

(2.6) ^(«),5 , "(Q)«-J2" , (0)«) = — 1 (r^ » A,(m))(F'(m), r t »r,(n)). 

2(Aj(u) - Xi(u)) 

If Aj is genuinely nonlinear, the left hand side of (2.6) is zero if and only if the rarefaction curve is a 
straight line, because r, • Ti(u) is orthogonal to r,i(u). 

The flux function / thus satisfies the following assumptions: 

HI) the eigenvalues Ai of Df are linearly degenerate or genuinely nonlinear; 
H2) the rarefaction curves form a system of coordinates; 
H3) shock and rarefaction curves coincide. 

The system ( [2 . 1| ) has thus nu linearly degenerate fields \i , corresponding to the Riemann invariants w. k , 
and n gn i = n — nid genuinely nonlinear fields Afe, corresponding to the Riemann invariants Wk- In the 
latter case we have • rk(u) — for all u G E. 

Remark 2.1. If fl C R 2 , then the rarefaction curves Ri(s)u always generate a system of Riemann coordi- 
nates. Thus our assumptions are satisfied by the following classes of systems: 

i) both eigenvalues are linearly degenerate; 

ii) one eigenvalue is linearly degenerate, the other genuinely nonlinear and the rarefaction curves of 
the latter are straight lines; 

iii) both eigenvalues are genuinely nonlinear and the system is of Temple class. 

The various situations are shown in fig. |l|. Case ii) corresponds to the traffic model considered in Q, 
while case i) corresponds to 2 x 2 chromatography. 

Given the two points u~ , u + G E, with coordinates u~ = u(w^ , . . . , w~) and u + = u(wf , . . . , u>+), 
with 7^ , consider the intermediate states u(u>i), where 

(2.7) u Q = w(u~), uj, = (w+, . . . wt +1 , . . . ,w~), i=l,...,n. 
For alH = 1, . . . , n, we denote with Vi(u~ , u + ) the vectors defined as 

(2.8) Vi(u~,u + ) — u{uji) - u(bJi-\), 




Figure 1. The various situations for a 2 x 2 system considered in Remark 2.1 



and we define Ti(u , u + ) as 

!v i {u~,u + ) U(b)i) - u(u)i-i) _ , 
= it 111- =?= W ■ 
\v i (u-,u+)\ |u(w<)-u(<*-i)| 4 _ 8 
ri(wi_i) = n(u)i) if tor = w+ 

where Ti{u) is the z-th eigenvector of DF{u). We assume that the vectors ri(u~ ,u + ) are linearly inde- 
pendent for all u~,u + G E. This condition is satisfied for data in a sufficiently small neighborhood of a 
given point u G H We denote also with {l l (u~ , «+), i = 1, . . . , n} the dual base. 

We now define an approximated semigroup of solutions S v on a set E v C £7. The construction is 
similar to the one in |Q. For any integer v € N, set 

(2.10) = {it G £ : Wi{u) G 2 _I/ Z,i = l,...,n|, 
and let D V,M be the domain defined as 

(2.11) D"' M = |n : R i — >E V :u piecewise constant and Tot.Var.(u) < m|. 

Given u G i? 1 ', we construct a solution by wave front tracking. We first define how to solve the 
Riemann problem [it~,u + ], with u~,u + G E v . 

The solution to the Riemann problem u~,u + is construc ted by piecing together the solutions to the 
simple Riemann problems where uji is defined in (2.7). If the i-th field is linearly degenerate, 

then [u>i—i, u>i\ is solved by a contact discontinuity travelling with speed \i(u)i). If the i-th field is genuinely 
nonlinear and wf < w~ , then [cj.;_i,u;i] is solved by a shock travelling with the Rankine-Hugoniot speed 
<Ti(uJi^i, LOi). Finally, if the i-th field is genuinely nonlinear and wf > w~ , then is solved by a 

rarefaction fan: if wf = + p- l 2~ 1 ' , pi G N, consider the states 

=Wi_i, = K,...,w i ti."' 1 r + ^2~ l/ ,w m ,...,w ? 7), 1= l,...,pi. 

The solution will consist of pi shock waves [wy-i, cjj,;], travelling with the corresponding shock speed 

At time i = we solve the initial Riemann problems of u. Note that the number of wave fronts is 
bounded by 2" •Tot.Var.(u). When two or more fronts interact, we solve again the Riemann problem they 
generate, and so on. It is easy to show that at each interaction at least one of the following alternatives 
holds: 

i) the number of waves decreases at least by 1; 

ii) the total variation of the solution u(t) decreases by 2 1 ~ ly , 
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iii) the interaction potential Q(t), defined as 
(2.12) Q(t)= ]T \a a \\ap\<M 2 , 

Q,/9 approaching 

decreases by 2~ v . We recall that two waves a a , o~p of the families k a , hp, located at points x a , Xp, 
are considered as approaching if x a < xp and k a > kp. 

This implies that there are at most a finite number of interactions, so that we can construct our approx- 
imate solution for all t > 0. Note that S^u — u(t) is a semigroup of solutions, but not entropic due to 
the presence of rarefaction fronts. 

If the i-th family is linearly degenerate, the i-th Riemann coordinate Wi(t,-) of the solution can be 
constructed by solving the semilinear system 



(2.13) 



(wi) t + Xi(u(t, x))(wi) x = 
Wi(0, x) = Wi{x) 



Since it is a piecewise constant solution, with a finite number of jumps, the broad solution to ( 2.13| ) is 
well defined §: if we denote with x(t, y) the solution to the ODE 

(2.14) x = \i(u(t,x)), x(0)=y, 
then the solution to ( 2.13] ) is given by 

(2.15) Wi(t,x(t,y)) =Wi(y). 

In the following sections we will consider the dependence on the initial data uq of the genuinely nonlinear 
Riemann coordinates Wk(t, ■) and the map h\(y) defined as 

(2.16) hl(y) = x t (t,y), 



where Xi(t,y) is the solution to (2.13). 



3. Estimates on the shift differential map 

In this section we prove some properties of the shift differential map. These properties are closely 
related to the structure of (2.1), i.e. the conservation form, the coinciding shock and rarefaction assump- 
tion, which prevents the creation of shock when two jumps of the same family collide, and the existence 
of Riemann invariants, which prevents the creation of shock when two jumps of different families interact. 

Consider a wave front solution u(t, •) of (2.1), and assume that the initial datum u(0, •) has a finite 
number N of jumps o~ a , located in y a : 

N 



u(0,x) = ^ a aX[y a ,+°o)(x)- 
a=l 

If £ a is the shift rate of the jump er Q , define u e (t, •) as the front tracking solution with initial datum 

N 

(3- 1 ) u 9 (0,x) = ^2o- a Xly a +0^,+oc)(x)- 

a=l 

In the following, we will use the integral shift function, defined by 

(3.2) v(t,x) = Bm|-iy u e {t,y)-u{t,y)dy\. 

If u{t, ■) has a shock erg, located in yp, and if £p is its shift rate, it is clear that the following relation 
holds: 

(3-3) a p£p = v(t,y p +) - v(t,y -). 

We first recall the following result in [[ll] , obtained using the conservation equation: 
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Figure 2. Interaction with a sheaf of shocks. 

Lemma 3.1. Consider a bounded, open region T in the t-x plane. Call o~ a , a = 1, . . . , N, the fronts 
entering T and let £ a be their shifts. Assume that the fronts leaving T, say a'p, (3 — I, ... , N' , are linearly 
independent. Then their shifts £'g are uniquely determined by the linear relation 

N' N 

(3.4) EW = E^- 

0=1 a=l 



Remark 3.2. As observed in Jll]], formula (3.4) implies that the shift rates of the outgoing fronts depend 



only on the shift rates of the incoming ones, and not on the order in which these wave-fronts interact 
inside T. In particular we can perform the following operations, without changing the shift rates of the 
outgoing fronts: 

01) switch the order of which three or more fronts interact; 

02) invert the order of two fronts at time 0, if they have zero shift rate. 

The second lemma is concerned with a configuration where a sequence of contact discontinuities inter- 
acts with a wave of another family. 

Lemma 3.3. Consider a family of parallel contact discontinuities a a , a = 1, ... ,7V of the i-th linearly 
degenerate family and a single wave-front a of the k-th family, k ^ i. Let ^ a and £ be their initial shifts, 
respectively, and let £' a> £' be their shifts after interaction. Assume that £ Q = £ for all a. Then after the 
interactions all the shift rates £' a of the i-th family have the same value £' and 

(o ^ b b _ £(A' — A) - g(A' - A) _ g(A'-A)-g(A'-A) 

where A, A and A', A' are the speeds of the shocks o~ a and a, before and after interaction, respectively. 

Proof. Define the vector v in the t-x plane as the shift of the first collision point. By a direct computation 
one finds 



(3.6) 



A-A' A-A 



Since all the incoming shock of the linearly degenerate family have the same speed A, by simple geometrical 
considerations it follows that the vector v is constant during all interactions (fig. ||) . Formula ( |3.5|) follows 
easily. □ 

Remark 3.4. Note that this lemma allows us to perform the following new operation, without changing 
the shift rates: 

03) replace a family of contact discontinuities o~ a of a linearly degenerate, all with the same shift rate 
£, by a single wave a = o a with shift rate £. 

In the next lemma we will show that the existence of Riemann coordinates w implies a strong relation 
among shocks of different families. 
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Figure 3. Vector relations among shocks. 



Lemma 3.5. Consider two adjacent jumps belonging to different families, Oi and o~j, i < j , located at 
Xi > Xj. Let a[, a'j be their strength after interaction. Then the following holds: 

(3.7) span{(7i, (jj} = span-jV^, a'j}. 

Proof. If £j, £j are the shift rates before interaction, and ^ after interaction, then (3.7) follows easily 
from the conservation relation 

(3.8) a& 



because by assumptions no waves of other families are generated. By condition (2.8) the conclusion 
follows. □ 



Remark 3.6. Note that the previous lemma implies that the conservation relation (3.£) is bidimensional, 
i.e. the shocks o~i, Oj and a[, a'j lie on a two dimensional plane (fig. |§|). We can obtain then an identity 



which relates the the strengths a with the speeds A: substituting (3.5) in (3.8), since £, £ are arbitrary, 
we get 

(3.9) a, (Aj -At) = a[ (Aj - A, ) + ^ (A'j - A, ) , 

- J (A i -A i ) = ^(Ai-A i )+^(A;-A i ). 

One can show that if a Riemann solver verifies ( ft .9] ) for all couple of waves i, j, then there exists a flux 
function / such that the wave front approximation is a weak solution to ( 2.1). 

An important property of the shift differential map for Temple class systems is the fact that a 
perturbation to the initial data, initially localized in [a, b], remains in the neighborhood of the set 
L)i[xi(t, a), Xi(t, b)], where Xi(t,y) is the solution of the i-th characteristic equation starting at y. We 
now extend this property to hyperbolic systems satisfying the hypotheses HI), H2), H3) of section ||. 

Consider N jumps cr Q , a = 1, . . . , N, of some linearly degenerate family i, located at x a and corre- 
sponding to the jumps c(a)ej in the Riemann coordinates w: 

(3.10) a a — u(w(x a —) + c(a)ei) — u(w(x a —)), 
for some constants c(a), a = 1, . . . , N. 

Definition 3.7. We say that the jumps a a defined in ( jOC| ) are in involution if 



(3.11) 



E 

a=l 



c(a) = 0, 



i.e. the initial and final Riemann coordinate Wi is the same: Wi(x\— ) — u>i(xN+). 

Note that, by the existence of Riemann coordinates, this relation does not depend on the positions 
and strength of the shocks of the other families. We can now extend Lemma 2 in fTl|| to our systems: 
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Figure 4. Computation of the shift rate. 



Lemma 3.8. Consider a wave front tracking solution u. Assume that there are N shocks o~ a 

i) either of the i-th linearly degenerate family in involution, 

ii) or of the k-th genuinely nonlinear family, 

and let x a (t), < t < T, be the position of the shock a a , a = l,...,N. Then it is possible to assign 
shift rates to all shocks such that £i = 1 and the shift of all fronts outside the strip T = {(t,x);t £ 
[0, T], X\(t) <x< xpf{t)} is zero. 

Proof. We consider only the case of linearly degenerate family i, since in the other case the proof is 
exactly the one given in . 

Let x a (t), a = 1, . . . , N, be the position of the shock o~ a of the i-th family in involution, and let Wi be 
the value of the Riemann coordinate at x\(t)— = xi(0) — . For w S E, define w as the projection of w on 
the hyperplane {wt = Wi\, and ii = u(w). 

We choose the shift rates such that 



(3.12) 



d 
~d9 



f u e dy= Z(i)(t)o-i(t) = c(t,x)(u(t, 

J — OO /A^ 



x) - u{t,x)), 



Xi(t)<X 



1, we need to prove that (3.12) 



where c(t,x) is a scalar function different from only in [x\(t), xpfft)], and we recall that u(t,x 
u(w(t, x)). 

By imposing the value £i = 1, i.e. c(0, xi(0)~- ) = 0, c(0,a;i(0)-F 
can be satisfied at time t — 0. We have two cases. 

1) If the jump o~i belongs to the i-th family and is inside [xo(0), xn(0)], then set £ = c(t, Xj —). 

2) If the jump m belong to the fc-th family with k ^ i, then by assumption (2.8) and by (3.7) there 
exists a unique shift and a unique constant c(0, x+) such that 

&o-i + c(0,x-)(u(0,x-) - u(0,x-)) = c(0,x+)(u(0,x+) - u(0,x)). 



Since we assume that the shocks are in involution, setting £jv = c(0, x n — ) we have that ( 3.12| ) holds at 
time t = 0: in fact the last jump has size u(Q, xjv(0)— ) — u(0, xn(0) — )■ 

We now show that this property is conserved for all t > 0. This follows easily from conservation 
and Lemma 3.5. The proof is exactly the same as in |pr| : we repeat it for completeness. Consider the 
interaction between two shocks o~i and Uj in the point (r, y), see fig. |Z[ By inductive assumption, we have 
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for the states ui, u m and ui that 

(3.13) ^2 (T ~t( T )£~<( T ) = c i( u i ~ «0i 

x-,{r)<y 

Ci(ui -ui)+<Ti£t = C m (u m - U rn ), 

Using conservation we have 

(3.14) + $4+ $4, 
so that for the new middle state we have 

(3.15) c t (ui - u z ) + a'^'j = c' m (u' m - u' m ) = c r (u r - u r ) - cr^,' 
and using Lemma 3.5 we conclude 

spanjw; - u^a'^'A P|span|u r - u r , cr-j = span{u m - u m ). 

The same r elatio n proves that they vanish outside L: in fact, assume for example that q = and j < i. 
Then from ( 3.15 ) we get 

^j^j c mi. u m ~ ^m)i 

which implies that c' m = 0. This concludes the proof. □ 

Remark 3.9. Note that for discontinuities of a linearly degenerate family all shift rates has the same sign. 
Note moreover that if no waves of other families are present, then we shift all jumps a a by unit rate 1. 
This corresponds to the case considered in Lemma |3.3| , i.e. to the substitution of a family of contact 
discontinuities with a single jump, whose strength in this case is by the involution assumption. 

Using conservation and the previous lemmas, we obtain explicitly the shift differential map at a give n 
time t. We recall that, given the states u~ ,u + e E, we denote with r.i (u~ , u + ) the vectors defined in (|2.9|) , 
and with l l (u~,u + ) its dual base. Let Pj(u~ ,u + ) be the projection operator on span{ri(u _ , w + ), i = 

i jY- 

j 

(3.16) Pj(u-,u + )v = ^2(l l (u',u + ),v)ri{u-,u + ), 

i=l 

where (•, •} denotes the the scalar product in W 1 . 

Given a point {t,x), with u(t,x) continuous in x, define x. t the intersection of the backward z-th 
characteristics starting at (t,x) with the real axis {(0,x)}, and for all (0, y) let j(y) the index such that 
x j(y) < 2/ < x j{y)-i> j{y) — I, ■ ■ ■ ,n+l, with xq = +oo and x n+ i = — oo. Without any loss of generality, 
we can assume that in (0, y) there is a jump a of the A:-th family. 

Define the points wi , w r £ E by 

!w(t,x) j(y) = 1 

(wi(0,y-), . . .,w J -(j,)_i(0,y-),w J - ( j / )(t,z), . . ,,w n (t,x)) 2 < j(y) < n 
w(Q,y-) j(y) = n+l 

(w(0,y+) j(y) = i 

w r (x,y) = < {w x {t,x), . . . ,w j(y) _ 1 (t,x),w j{y) (0,y+), . . .,w n {0,y+)) 2 < j(y) < n 

[w(t,x) j(y)=n+l 
Moreover define the point w m € E by 
(3.18) 

w m (x,y) = (wi(t,x), . . .,w k (Q,y+), . . . , x), w j{y) (0, y+), . . .,w n (0,y+)) 2 < j(y) < n, 

if k < j(y), and in a similar way, if k > j(y), 
(3.19) 

w m (x,y) = (wi(t,x), . . . ,w j (y)- 1 (t,x),w j ( y )(0,y+), . . .,w k (0,y-), . . . ,w n (0,y+)) 2 < j(y) < n. 
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Figure 5. Wave pattern for the computation of formula 3.16 



Define P(x, y) as the vector 
(3.20) 

f o j(y) = i 

P(x,y) = I Pj( y )-i(wi,w m )a + Pj( y )-i(w m ,w r )(a - Pj^_i(wi,w m )a) 2 < j(y) <n + l,k< j(y) 
[ P j(y)-i( w h w m)(Pj(y)-i(w m ,w r )a) 2<j(y) <n,k>j(y) 

where wi = wi(x,y), w m — w m (x,y), w r — w r (x,y) and a is the initial jump in (0,y). Consider now a 
front tracking solution u & , obtained by shifting the initial jumps <j a in y a with rates £ Q . 



Theorem 3.10. Ifv(t,x) is the integral shift function ofu (t,-), defined in (3.2), then 

1 



(3.21) 



v(t, x) = lim < — - 



u (t,y) - u(t,y)dy 



Proof. The theorem will be proved outside the times of interaction, because the Lipschitz dependence in 



L 1 of the approximate semigroup implies the validity of (3.21) for alH > 0. 

If is sufficient to show that P(x,y a )t; a is piecewise constant, with jumps only at the points xp 
where u(t, •) has a shock crp, and the following relation holds: 



(3.22) 



E 



(P(xp+,y a ) - P(xp-,y a ))£ a = o-ptp, 



lim 



P(x,y a )£ a = 0, 



where is the shift rate of 073, located in x. Note that by ( 3.2C ) the second equality of ( 3.22 ) is trivially 
satisfied. 

By linearity in the shift rates £ Q , we can consider the case in which a single shock is shifted, let us say 
a at y: (3.21) becomes 



(3.23) v(t,x)=P(x,y)t 

Formula ( 3.16 ) follows from the following considerations: consider a wave front pattern, fig. ||, where for 



simplicity we assume that k < j(y). The states wi, w m are computed considering the Riemann problem 
generated by adding to the fc-jump a in (0, y) all the i-waves starting from the left of (0, y) and ending 
in the right of (t, x) and all the ?-waves, with i ^ k, starting from the right of (0, y) and ending in the 
left of (t, x). The jump w m ,w r is a single wave of the fc-th f amil y formed by adding all the k- waves 
bet ween (0, y) and (t,x). Using the definition of v(t, x) given (^), one obtains easily the second case 
of ( 3.20 ): in fact the shift rates of the shocks in the left of (t, x) is given by the shift rates of the jumps 
of the Riemann problem wi,w m ending in the left of (i, x), Pjf y -\—i(wi,w m )cr, plus the shift rate of the 
shock w m ,w r , Pk(w m ,W r )(a — Pjt y -\_i(wi,w m )a). Since only the i-waves with % >j(y) > k are present in 
<t — Pj( y )-l(wi,w m )a, then P k (w m , w r )(a - P j(y ^ 1 (wi,w m )a) = Pj( y )-i(w m ,w r )(cr - P j ^_ 1 (wi,w m )a). 
The other cases can be computed in a similar way: in this case one solves the Riemann problem w m , w r 
in (0, y), and consider the fc-wave Wi, w m starting in the left of (0, y) and ending in the right of (t, x). 
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Figure 6. Computation of the shift rate in the case of Lemma |3.11 



From the above considerations it is clear that P(x,y) is piecewise constant, with jumps only when in 
(t, x) there is a i-shock a': in fact otherwise the wave front pattern used to compute P(x, y) remains the 
same. Let {z p : p = 1, . . . , M} be the set of the starting points of all shocks arriving in (t, x), and define 

(3.24) z~ " 



mm z. 
p 



We consider two cases: 

1) the shocks arriving in (t, x) start on both sides of (0,y): z~ < y < z + . In this case, (P(x+,y) — 
P(x—, y))t; is the shift rate of the i-shock starting in the Riemann problem wi(x—, y), w m (x+, y) if 
i > k (w m (x—,y),w r (x+,y) if i < k) which collides with a fc-shock w m (x+, y), w r (x+, y) (wi,w m 
if i < k): in fact the only difference is that in w m (x—,y),w m (x+,y) there is a shock of the i-th 
family starting in (0, y), and i is genuinely nonlinear. Finally, using • ri{u) = and Lemma 3.5, 
one can change position to the i-wave and the remaining /c-wave w m , w r , whose strength does not 
change. 

If i = k, there are no fc-shocks starting on the right (left) of (0, y) and ending on the right (left) 
of (t, x), so that (P(x+,y) — P(x—,y))(, is the shift rate of the i-shock of the Riemann problem 
wi(x-,y),w r (x+,y). 

2) the shocks of the i-th family arriving in (t, x) start either in (— oo,y) or (y, -t-oo): assume for 
definiteness that y < z~ . In this case the difference (P(x+,y) — P(x—,y))£ is the shift rate of the 
shock a' colliding with the shifted shocks of the Riemann problem wi(x—,y),w m (x—,y) in (0, y), 
crossing the jump w m (x—, y), w r (x—, y), and finally overtaking a' . In fact one can use Lemma p.5| 
(and ri • ri(u) = if i is genuinely nonlinear) to obtain the wave pattern of fig. |^. 

The various cases will be proved in the following lemmas. 



Lemma 3.11. 

e is 

(3.25) 



Assume that z <y< z + , i.e. case 1). If the shock a' is of the i-th family, then its shift 
£V= (P(x+,y)-P(x-,y))Z. 



Proof. We follow closely the method of |TT[. Assume for definiteness k < j(y), the other cases being 
similar. The basic idea is to reduce the computation to the single Riemann problem wi(x—, y), w m (x+, y), 
with eventually a single fc-wave w m ,w r . 

Consider fig. [| By Lemma [O], we can simplify the wave configuration considering only the fronts 
crossing starting in the right of (0, y) and ending in the left of (t, x): in fact we can move the other fronts 
to ±oo without changing the shift rate of a'. 

We can now shift the initial position of the waves of the i-th family merging in x such that their initial 
position coincide with y, without changing the shift rate This operation can be repeated for all shocks 
of genuinely nonlinear families. 

Finally, we can move the shocks of the linearly degenerate families such that they have the same 
sequence of interaction with the other shocks. This means that, if x\ is the position of the j-th shock of 
the i-th linearly degenerate family, the only interactions among shocks occurring in the sector \x\(t\ xf(t)] 
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Figure 7. Computation of the shift rate in the case of Lemma |3.12 



are the one involving one i-th. wave and one fc-th wave, with k i. Using Lemma 3.3, we can at this 
point substitute them with a single shock, whose strength is the sum of the strengths of the i-waves. 
Finally we move their position at t = such that it coincides with y: we obtain the wave patterns of 
fig. |^. To conclude, we just need to prove that the Riemann problem obtained in this way is exactly 
Wi(x— , y),w m (x+, y) and that the remaining £;-wave is w m (x+, y), us r (x+, y). 

By the previous argument, the strength of the shock of the j-th family j < i, j ^ k, is given by the 
j-waves starting in the right of y and ending in the left of x: since they are the only j-wave crossing the 
segment [{0,y— ), (t,x+)], it follows 

wi,j(x,y) = wj(0,y-), w r> j(x,y) = Wj(t,x+). 

The other relations for j = k and j > i follows in the same way. Finally, for j = i the jump is 
Wi(t,x+) — Wi(t,x—). Note that the wave pattern is the same obtained in 1). □ 

We consider only the case y < z~ , since the other is entirely similar. 

Lemma 3.12. Assume that y < z~ . Then the shift £' of a' is given by 

(3.26) ?°'=(P(x+,v)-P(x-,v))Z. 

Proof. The hypothesis implies that the i-th shocks ending at x starts in the right of y. With the same 
simplification considered in Lemma 3.11, we reduce to the Riemann problem Wi(x, y), w r (x, y) in y, such 
that the waves of the j-th families, j > i, generated at y collide with the i-wave in x a (see fig. Q), 
after overtaking the fc-wave w m ,w r . The conclusion follows easily, since the wave pattern is the same 
considered in 2). □ 



This concludes the proof of Theorem 3.1C 



□ 



Finally we extend to our case the following result proved in ]il[ 



Proposition 3.13. Let u be a wave-front tracking solution, and consider two wave-fronts, x(t) and y(t), 
t G [0, T] . Then there exists a second front tracking solution u such that the initial and final positions of 
the two shocks is the same, and Tot.Var.(it) is uniformly bounded. 

Proof. For genuinely nonlinear fields, the proof is the same as in pT| . We then restrict the proof to the 
case of a linearly degenerate fields i. 

Assume that there exists two jumps o~\, o~i of the z-th family, with positions z\(t) < z-x (t), such that 

(3.27) x(0) £ MO), z 2 (0)] and y(0) £ [^(0), z 2 (0)], x(T) £ [ Zl (T), z 2 (T)} and y(T) £ [ Zl (T), z 2 (T)}. 
For definiteness, assume u>,(0, Zi—) < u>i(0, and the following conditions is satisfied: 

(3.28) u><(0,*i-) G [w l (Q,z 2 -),w(0,z 2 +)}. 

Let cr Q , a = 1, . . . , N be the jumps of linearly degenerate family i in the strip [zi(0), z 2 (0)): if we define 

cttv+i = u(wi(0, zi-)) - u(wi(0, z 2 -)), 
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Figure 8. Cancellation among contact discontinuities. 



it is easy to verify that the shocks er a , a=l,...,JV+l, are in involution. By Lemma 3.8, we can then 
moving the jumps to the left until either Zi(t) meets the wave fronts x(t), or z\(t) coincides with another 
shock of the i-th family (fig. ||). It is clear that we can repeat the same procedure also in the following 
cases: 

i) Wi(0, zi—) > Wi(0, Zi+) and Wi(0, z\—) E [u)<(0, z 2 +), w(0, z%— )]; 

ii) Wi(0, z%— ) < Wi(0, z 2 +) and Wi(0, z 2 +) G [wi(0,zi+),w(0,zi—)]; 
hi) Wi(0, Z2— ) > Wi(0, z 2 +) and Wi(0,z 2 +) <E [wi{0,zi—),w(0,zi+)}. 

It is now easy to prove that the total variation of the jumps of the i-th family satisfying ( |3.27 ) can at 
most be SHioH^,. Since x(t), y(t) divide the lines t = and t = r in three regions, the total variation of 
u>i is bounded by 27||w;|| 00 . □ 



4. Estimates on characteristics 

In this section we prove some estimates on the solution Xi(t, y) of the characteristic equation: 

f .,, j ii = Xi(u(t,Xi)) 

[ ' \ Xi(0) = y 

We assume for simplicity that the i-th family is linearly degenerate, however the same results are valid 
for characteristics of a genuinely nonlinear family if the following condition holds: for all r there exists 
an e such that in the strip {(t, x);t < t < T,Xi(t,y) — e < x < Xi(t,y) + e} there are no shock waves 
of the i-th family. Given front tracking approximation u, Xi(t,y) is unique, since it crosses only a finite 
number of transversal jumps, and it depends Lipschitz continuously on the initial data y (see [£)). 

We want to give uniform estimates on this dependence. The idea is to suppose that in y there is a 
shock cr e of the i-family of size e: Wi(0, y+) —Wi(0, y—) = e. Since by assumption no shocks of the i-family 
collide with cr e , it is easy to construct a wave front solution: for x < x(t,y), the solution it c (t, •) takes 
values in 

E u '~ = I?/ : Wj(u) £ [dj, bj] fl 2~ ly Z, j = 1, . . . , n\, 

while for x > x(t,y), enlarging E and assuming e sufficiently small, 

E u ' + = iu : w 3 {u) G [aj,bj] n2^Z,j ^ i,w t (u) G [a t ,b t ] n {2 _1, Z + e}|. 

The following lemma proves the continuous dependence of the solution u e (t) and the position x\(t, y) 
of the shock a e w.r.t. e. 

Lemma 4.1. Consider a fron t tr acking solution u, with initial data uq and the characteristic lines 
Xi(t, yi) < Xi(t, y 2 ), defined in (4.1) for a linearly degenerate family i. Let u e the wave front solution with 
initial data u(wq), where Wq is defined as 

!w(u (x)) x < yi 

w(uq(x)) + ea yi < x < y 2 
w(uq(x)) x > y 2 
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Then there exists constants L, L 1 , depending only on the total variation of the initial data uq, such that 
for allt>0 



(4.3) 



\u(t, x) — u e (t, x)\dx < Le 



where x\(t,yj) is the position of the shock crj starting in (0,yj). 

Proof. The first inequality is an easy consequence of the L 1 continuous dependence for front tracking 
solutions, see §]. For the second one, note that all the shocks different from a e have size uniformly bigger 
than 0, so that their position is shifted of the order e. Thus the second inequality follows by standard 
ODE perturbation estimates, see [jj. □ 



An easy application of the previous lemma together with Proposition 3.13 implies that to compute x\ (t, yi) 
and 22 (i, 2/2), we can ac tual ly consider in (4.1) a solution u with uniformly bounded total variation, so 
that the constant L' in ( |4.3| ) is independent on the total variation of Uq. 
We now estimate the dependence of Xi(t, y) w.r.t. u. 

Proposition 4.2. Let £ Q be the shift rate of the jump o~ a in u(0, ■), and denote with xf the solution to 

" i? = Ai(u fl (t,^)) 

*?(o) = y 

where u e (t) is the shifted front tracking solution. Then there exists a constant D independent of the total 
variation of u such that 

( 44 ) i:„ x e i {t,y)-x i {t,y) 



lim 

6>^0 



Proof. We give a sketch of the proof, for details one can see |Tl| . 

If e is the size of the shock a e located in (0,y), then we can apply Theorem 3.1C to compute its shift 
£ e : by formula (3.21) we obtain 



(4.5) 



f a e = ^2(P(x+, y a ) - P{x-,y a )% 



If 9 is sufficiently small, then we have 



Xi' e (t, y)-xf(t, y) 



where x i ' e (t,y) is the position of the shifted shock and x1(t,y) is its original position. Note that 
(P(x+, y a ) — P(x—, y a ))£,a is the shift rate of the shock <r e , after colliding with the shocks of the Riemann 
problems wi,w m and w rn ,w r . Their total shift is proportional to |<7q,£q,|, and after the interaction with 
er e , the shift of the latter is proportional to | cr e 1 1 ct^^q, | . Thus taking the limit as e tends to of (fi~5|), we 
obtain for e sufficiently small 

x 6 i{t,y) ~Xi{t,y) 



which implies (4.4). 

We prove now the uniform Lipschitz continuity of the map y 



□ 



Xi(t, y) for alH > 0. 



Proposition 4.3. Consider two characteristic lines Xi(t,yi) Xi(t,y2), solution to (4.1). There exists 
C > 0, depending only on the system and the set E, such that 



(4.6) 



1 < xj^m) -x}(t, yi) < c 
C ~ y 2 - yi 



Proof. As in the previous proposition, let e be the size of the shock cr e (t) located in (0, y) in Riemann 
coordinates. If £(i) is its shift rate, then for 9 sufficiently small by Theorem 3.10 we obtain 

xf(t,y + 9Z)-xi(t,y) 



(4.7) 



-a'{t) =C (t)o- e (t) =r i (w l ,w r )(l\w l ,w r ),o-(0))Z(0). 
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Figure 9. Decay of positive waves. 

In fact, by assumption, in the simplified wave patterns to compute the shift rate of a e , there are no waves 
of the i-th family different from o~ e . Dividing by e and taking the limit as e tends to 0, we obtain 

Xi(t,y + 6£) - Xi(t,y) d d ; i , 
6 dw ■ ' = duy ^ ^ (wi,w r ), n(0, y))£{0), 

which implies 

We use the fact that a e (t)/e tends to du(0, y)/dwi-ri(wi,w r ) as e — > 0. Since E is compact, the conclusion 



(4.6) follows easily. □ 



Remark 4.4. The above proposition implies that the map h\ defined in ( [2.16 ) is uniformly Lipschitz, 



independent on the total variation of uq, together with its inverse map (/i*) 1 



To end this section, we give a different proof of the following result given in 1 1 



Proposition 4.5. If x{t), y(t) are the positions of two adjacent k-rarefaction waves, then for some 
constant k > one has 

(4.9) y(r) - x(t) > ktT v . 

Thus for all r > the total variation of the Riemann invariant u>k of the k-th genuinely nonlinear family 
with N shocks at t — is bounded by 

(4.10) Tot.Var.{ Wfc (r, ■); [a, b}} < 2(6 ~ a) + \\w k \\ Loo + (N + 1)2 1 -". 

Proof. Consider two adjacent A;-rarefaction fronts x{t) and y(t), and let t a , a = 1,...,N, be the in- 
teraction times of x(t), y(t) with other waves in the interval [0, r]. Fixed ti G (ta,t a +i) for some a, 
let z(t, x(ti)) be the characteristic line of the fc-th genuinely nonlinear family starting in (ti,x(ti)) (see 
fig. H). Assume t i+ i > ti sufficiently close to t% such that S (ia,*a+i) and z(t, x(U)) does not collide 
with shocks of other families for t G [tj, t^+i] . Let z(t, a;(tj+i)) be the characteristic curve starting in 
(ti+i, x(ti+\)). By the assumption of genuinely nonlinearity, at time ij+i we have 

z(t l+1 ,x(U)) - z(t i+1 ,x(t i+1 )) > c(t i+ i -t^"^ 1 , 



for some constant c, depending only on E. Using Proposition 4.3, at time r we have 



(4.11) z(t,x(U))- z{ T ,x(t l+l )) > ^(t i+1 -t i )2-' / - 1 . 

Repeating the process, it is possible to find a countable number of times ti such that 

lim ti=ta, lim ti = t a +i, 

i — > — oo i — >+oo 
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and using ( 4.11 ) we get 

(4.12) z(T,x(t a )- Z {T,x(t a+1 )) > ^(tc t+ i-t & )2- / - 1 . 



Repeating the process for y{t) and for all intervals {t a+ i,t a ), we obtain ( |4.9| ) where k — c/C. 

The second equation follows noticing that the total amount of positive jumps in the interval [a, b] is 
bounded by (1 + N)2~ u + (b - a) /kt. □ 

5. Proof of the main theorem 

In this section we construct the semigroup S on L°°(R;E). In (2| it is shown that for all M, there 
exists a semigroup S M defined on the domain 

(5.1) D M = ju : R >-> E : Tot.Var.(u) < Afj, 

which is the only limit of the wave front tracking approximations constructed in section |^. We study now 
the dependence of the solution on the initial data u € D M . We consider separately the case for genuinely 
nonlinear and linearly degenerate families. 

Proposition 5.1. Consider a front tracking solution u, such that u(0, •) has N jumps a a , a = 1, . . . ,N, 

and let £ Q be their shift rates. Given r > 0, denote with ap the jumps in the Riemann invariant u>fc(r, •) 
of the k-th genuinely nonlinear family. Then there exists a constant K , depending only on the system 
and the domain E such that 

N 

(5-2) < ^(1 + ^2-") Y,\^o\- 

P a=l 

Proof. The proof follows by Theorem 3.10| and Proposition 4.5. In fact, fixed a shock a s , using Theorem 
3.10, we have that at time r for a shock 17,3 of the i-th family there exist D 1 

(5.3) l^o-pl < D'\£ a aa\, 



if the shock ap starts on both sides of cr a , or, using the same estimate of Proposition 4.2, 
(5-4) \U<D\a & i & l 



if ap start on one side of cr a . Since there is at most 1 shocks such that (5.3) holds, and the interval of 



influence ] , using Proposition |4.5| together with (|5.3|) and (pA) we obtain 



E 



Zpc-p\ < D'\£ & o- & \ + D\a a ^\ ■ Tot.Var.{ Wfc , [x 5 - \t,x & + At]} < F(l + 2-^ a a 6 



The conclusion follows the linearity of the shift differential map. □ 
Using the results of the previous section, the following result is trivial: 

Proposition 5.2. Consider a wave front solution u, such that it(0, ■) has N jumps a a , ct — 1, . . . , N, 
and let £ Q be their shifts. Consider the equation ( |4.l| ) ; whit the eigenvalue Xi linearly degenerate. Fixed 
t > 0, then the shift & of Xi(r,y) is bounded by 

N 

(5.5) \^\<DY / \^<J Q \- 

a=l 



Proof. This is a corollary of Proposition 4.2. □ 

Using the above propositions, we can prove the following theorem: 

Theorem 5.3. Consider two initial data u\ and U2, and denote with uij.k{t,-), the k-th Riemann co- 
ordinate of S M Uj, j — 1,2, corresponding to the k-th genuinely non linear family. Moreover, let hZ 



j = 1, 2, the map defined in (4.1) for the i-th linearly degenerate family. Then there exists a constant K 1 , 
independent of M , such that the following estimates hold: 

(5.6) / I wi : k(t, x) — W2.k (t, x) \dx < K 1 / \ui(x) — U2(x)\dx, 
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(5.7) 



sup | h\ j (x) — h\ i (x) | < K' I \ui(x) — U2(x)\dx. 

t>0,x£E 



Proof. Consider two piecewise constant initial data u\, u 2 in D M ' V , and construct a pseudo polygonal 
path 70 : 9 \ — > Ug, connecting u± and u 2 , such that 



||7o|| L i < E\\u% - v%\\ L1 . 

We can assume that Ug has a finite number N of jumps. If we denote with 7^ the path 9 
have by Proposition pTft 



S"v,g, we 



(5.8) 



\w2,k( T ) -Wiu(t)\ 



L 1 



< 



{Y T ) k ri < ^(1+^2^)11701 



L 1 



<if'(l + A^2" 1/ )|U 2 -iii| 



If now f — > +00, since k (r) converges to u^a^t), we obtain (5J:). Since this estimate does not depend 
on the number of initial jumps N, we can extend it uniformly on D M . 

Using the same pseudo polygonal path, in a similar way we can prove that 

\ x 2,ii T >v) ~ ^i,i( T >y)| < K'\\u 2 - Ui\\ L1 . 

This shows that x"(t, •) converges uniformly to the solution Xi(r, •) as v 
implies that 

\x2,i(r,y) - x lti (T,y)\ < K'\\u 2 

This concludes the proof. 

We can now define S on the domain L°° (R; E) : 
Definition 5.4. For all u e L°°(R, E), let u M G D M be such that 
(5.9) 



L 1 ' 



-00 and u" — ► u. It also 



□ 



lim u M 

M— >+oo 



tt in Ljk c . 



Define <S t u as 
(5.10) 

where the limit is in Ll 



S t u 



lim S t M u, 



It is easy to prove that the right hand side of ( 5.1C ) is a Cauchy sequence in every compact set [a, b): 
in fact, using the finite speed of propagation, we can consider u with compact support [a — At, b + At]. 
For the components Wk of the &-th genuinely nonlinear family, it follows directly from (5.6), while for a 
linearly degenerate component Wi, let w be a Lipschitz continuous function such that 

u>j(0, x) — w(x) \dx < e. 



By Theorem p . 3| we have for u\, u 2 6 D such that \\u — u^l 1 < 8, i = 1,2 



sup |/i*i,i(x)-/iL(ic)| <!£■'<$, 

t>0,xGR 

and it follows by easy computations that 

(5.11) \\wi <i (t)-W2,i(t)\\ L1 <\\wi li (t)-wo(h\ ii ) \\ L1 + ||ty 2 ,i(t) -<5o (/4,i) |Li + 

< C||w M (0) - w|| Ll + C||«;2,t(0) - w|| Ll + L(b - a)G\\u 2 - 

< 2C(e + S) + L(b-a)GS, 

where L is the Lipschitz constant of ui. This shows that wf 4 (t) is a Cauchy sequence for all t > 0, because 
the right hand side of (5.11) can be made arbitrarily small. We can now prove the main theorem: 



"1 



Theorem 5.5. The semigroup S : [0, +00) <g> L°°(R;i;) 1 — ► L°°(R;E) defined m (|5.10|) is the only 
continuous semigroup on L°°(R;E) such that the following properties are satisfied: 
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i) for all u n , u G L°°(R; E), t n , t G [0, +00), with u n — > u in L} oc , \t — t n \ — » as n — > +oo, 
(5.12) lim S tn u n =S t u in L} oc ; 



n — »+oo 



ii) each trajectory t t— * iSt^o is a weafc entropic solution to the Cauchy problem 

u t + f(u) x = 



(5.13) 



it(0, x) = uq(x) 



with u G L°°(R;E); 

hi) if uq is piecewise constant, then, fort sufficiently small, Stuo coincides with the function obtained 
by piecing together the solutions of the corresponding Riemann problems. 

Proof. The statement follows easily, since we proved that S t u is the unique limit of wave front approxi- 
mations, and for data with bounded total variation we can apply the results in □ 

Remark 5.6. Note that what we also proved that the characteristic equation ([O]) is well posed for L°° 
data: the solution Xi(t,y) is Lipschitz continuous w.r.t. both variables. This is not trivial, since even for 
2x2 systems not in conservation form the dependence is Holder continuous, while for general n x n the 
solution does not exists. 

Note moreover that semigroup S is continuous, but not uniformly continuous. However if the initial 
data takes values is a compact set of L 1 n L°°, then the semigroup becomes uniformly continuous. This 
extend the Lipschitz continuity when the initial data have bounded total variation. 
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